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Abstract. Observational data show that the correlation between supermassive black holes 
(Mbh) and galaxy bulge (.Mbulge) masses follows a nearly linear trend, and that the correlation 
is strongest with the bulge rather than the total stellar mass (_M ga i). With increasing redshift, 
the ratio Y — A^bh/A 7 ! bulge relative to z = also seems to be larger for A4bh > 10 8 ' 5 Mq. This 
study looks more closely at statistics to better understand the creation and observations of the 
Mbh — Albuigo correlation. It is possible to show that if galaxy merging statistics can drive the 
correlation, minor mergers are responsible for causing a convergence to linearity most evident at 
high masses, whereas major mergers have a central limit convergence that more strongly reduces 
the scatter. This statistical reasoning is agnostic about galaxy morphology. Therefore, combining 
statistical prediction (more major mergers =>■ tighter correlation) with observations (bulges = 
tightest correlation), would lead one to conclude that more major mergers (throughout an entire 
merger tree, not just the primary branch) give rise to more prominent bulges. Lastly, with re- 
gard to controversial findings that F increases with redshift, this study shows why the luminosity 
function (LF) bias argument, taken correctly at face value, actually strengthens, rather than 
weakens, the findings. However, correcting for LF bias is unwarranted because the BH mass 
scale for quasars is bootstrapped to the A4bh — o"* correlation in normal galaxies at z = 0, and 
quasar-quasar comparisons are mostly internally consistent. In Monte-Carlo simulations, high T 
galaxies are indeed present: they are statistical outliers (i.e. "under-merged") that take longer 
to converge to linearity via minor mergers. Another evidence that the galaxies are undermassive 
at z > 2 for their Mbh is that the quasar hosts are very compact for their expected mass. 
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1. Introduction 

T he discoveries of fundamental correlations betwee n A^bh with stellar velocity disper- 
sion ( Ferrarese fe Merritd2000t Gebhardt et al.|[2000h and M bulge ( Kormendv fe Richstone 
19951: iMagorrian et al.lll998l) have led to modern views th at supermassive black hole ac- 



tivities may have a strong impact on galaxy evolution (e.g Kauffmann fe HaehneltlliOQfj 
Granato et al.ll2004lDi Matteo et al.ll2005l) . The correlation between Mbh and A^buigc 



is remarkable in that it is almost linear, has only a scatter of 0.3 dex, and holds true over 
5 orders of magnitude in A^bh dynamic ran ge. Locally, the ratio of M. bulge /A^bh ~ 800 
(|Marconi fc Hundl2003UHaring fe Rbd^OO^ . 

How did the correlations come about and how do selection biases affect our obser- 
vations of the correlations? Direct cause and effect are not only possible, there are nu- 
merous theoretical propos als. While quasar feedback is one of the most investigated 
and fa vored explanations ( Di Matteo et al. 2005 : Robertson et al.l 20061 Hopkins et al 



2007a), galaxy mergers may perhaps share a role. This study therefore isolates the role 
of merger statistics to examine how it might affect the growth of the A^bh — Albuigc cor- 
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Figure 1. How mergers cause Mbh and M ga i to correlate, a) No correlation initially. Objects 
at location 1 increases .M ga i more quickly than Mbh upon merging, compared to objects at 
location 2. The correlation steepens over time due to minor mergers, b) Steep correlation initially. 
The opposite situation occurs, i.e. .M ga i grows more quickly than .Mbh for objects at 1 compared 
to 2. Symmetry between scenarios a and 6 produces linear, relation, asymptotically. Panels on 
the right are Monte-Carlo simulations with the initial conditions shown on left. 

relation. It also examines more closely how luminosity selection biases affect the inference 
of .Mbh — -Mbuigc correlation and its evolution since z s=s 4, as inferred using quasars. 



2. Merging Statistics: How .Mbh — .M bulge Correlation Can Result and 
Why it is Important to Also Consider the .Mbh — -Mtotai Relation 
How galaxy mer ging affects the BH vs. bulge correlations h as been considered by 
several studies (e.g. Ilslam et al]|2003t ICiotti fc van Albadall200l from a purely statisti- 
cal standpoint, and using specific initial conditi ons and ass umptions (e.g. no scatter or 
starting with a prior correlation). Going further, [Peng] ( 2007 ) shows the two most salient 
features of the Mbh — .Mbuige correlation: linearity and strong correlation with bulges, 
can be attained without having t o mak e assu mptions about the initial conditions. The 
heuristic toy model proposed by IPend (l2007lh shown in Fig. [T] explains that a linear 
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Figure 2. Effects of major mergers based on Monte-Carlo simulations. The solid line corre- 
sponds to a linear correlation; it is not a fit. The contours show the distribution of points at 
the end of the Monte Carlo simulation, a) No major mergers iV ma j = 0. The dashed line shows 
roughly the mean of the points with slope fixed to a linear correlation, b) Cumulative major 
mergers < iV ma j < 10 over all progenitors in an merger tree of a galaxy, c) Cumulative major 
mergers -/V ma j ^ 10 over all progenitors galaxies. 



quality of the A^bh — A'Ibuigc relation comes about because of minor mergers. Minor 
mergers affect the A^bh _ Af ga i relation in a way that over some number of events would 
attain a cosmic average ratio in A'fBH'/Algai. That there is necessarily a correlation can 
be reasoned from symmetry argument; that the correlation trends toward linearity can 
be understood by noticing that the only way minor mergers can no longer change the 
-AfBH/Algai ratio is when it is the same value everywhere along the mass sequence. 

Major mergers, however, play a different role: they do not significantly change the ratio 
of Msa/Mgai because the galaxies involved in merging have both sim ilar {Mbh) and 
(A'fgai) by definition of major mergers. As explained in more detail in IPenel (j2007h . m 
the limit where major mergers are occurring between identical mass galaxies, the scatter 

uIAAbk i+A^bh 2) 

of log(AI B H) after merging roughly goes like a (log (/x merge )) = ^bh i+A1bh 2) ' due 
to the central limit theorem. Thus, because the sum, (A^bh.i + A!bh,2)j increases more 
quickly than the dispersion, a (A^bh.i + A4bh.2)i the scatter in the AIbh distribution 
after merging, a (log (/i mer ge)), decreases with each major merger. The effect of central 
limit convergence due to major mergers can be quite dramatic. Fig. [5] shows results from 
one possible Monte-Carlo simulation, in which there is no correlation initially, and has 2 
orders of magnitude in scatter. Minor mergers (Fig. ^jpf alone do not reduce the scatter 
significantly by the end of the simulation. However, major mergers cause a rapid decrease 
in scatter in only a few events. Note that the relevant accounting of major mergers is the 
cumulative sum over the entire merger tree, i.e. over all progenitors, their progenitors, 
etc., rather than the more common approach of tracking the main branch. 

It is worth noting that statistical reasoning docs not predict morphology from first 
principles. Therefore statistical reasoning can explain the observations of a tight A^bh — 
M bulge correlation if and only if massive bulges were preferentially formed through more 
major mergers than disky galaxies, summed over all progenitor histories. Even though 
the notion that major mergers lead to formation of bulges is now widely regarded to be 
true, it is interesting that it can be reasoned purely from statistical principles and the 
known existence of a tight A^bh — A'fbuigc correlation. Furthermore, the A^bh — A^buigc 
relation may be a special case of the AIbh — Aftotai relation, despite the latter having a 
larger scatter. Thus, to understand the co-evolution of galaxies with A^bh, one ought to 
consider both the A(bh — Af bulge and AIbh — Aftotai relations. 
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Figure 3. Observed z > 1.7 quasar data (open circles) from iPeng et al.l (|2006al ). Solid circles 
are z = normal galaxies. The absolute luminosities are how the galaxies would appear after 
accounting for luminosity fading. 

Figure 4. The intrinsic correlation. What it means for the intrinsic correlation between 
Mbh and .Mbuigc to be linear. This distribution is also called "the prior" and the condi- 
tional P(A4BH|A4bui ge ). To simplify discussion, we assume P(A / lBH|A4bui gG )=P(A1buige|AlBH) 
as shown; doing so does not affect the conclusion qualitatively. 



3. Luminosity Function Bias of Galaxies, Quasars and Other Biases 

In recent years, there have been a number of efforts to study the .Mbh — M. bulge 
correlation beyond the local universe u si ng quasars, radio galaxies, and oth er means (e.g. 
McLure et al.ll2006t IPeng et al.l[2006aM b1: 1 Woo et al-bood llreu et alJl2007f ). For the most 



part, the studies find that the central BHs were larger at z > 2 in the past for a given 
bulge stellar mass by a factor r of 3 < T < 6 (|Peng et al.ll2006alfbh . shown in Fig.H These 
findings have been called into question by other studies on the basis that the luminosity 



function (LF) selection was not explicitly accounted (e.g. Lauer et al. 20071 ). However, 



the criticisms have not been very germane, both because the BH mass scale in quasars is 
not on an absolute scale, and the LF bias goes in the opposite direction in quasars than 
claimed, as discussed below. The issues are subtle and have led to substantial confusion. 



3.1. Revisiting the Luminosity Function Selection Bias to See Why It Affects the 
■Mbh — .Mbuigc Correlation in Galaxies Differently From Quasars. 
Figures 2HS] illustrate schematically the idea of the LF selection bias. Figure 2] is the 
prior that there is an intrinsic, perfectly linear, correlation between .Mbh and .Mbuigc- 
This intrinsic correlation P(.MBH[.Mbuige) is also known as the conditional. Even though 
a linear relation does not require P(.MBH|-Mbuigc) to be the same as P(-Mbuigel-MBH), 
doing so in the discussion below does not affect the directional se nse of the conclusions. 



The luminosity function bias was pointe d out a t least as early as I Adelberger &: Steidel 



(|2006l) ; ISalviander et all (|2006t ; lLauer et a" 



(|2005|h and more recently by iFine et al 
(|2007l ). In essence, the act of selecting a sample of galaxies to observe leaves an imprint 
of the LF on the correlation of .Mbh vs. .Mbuigc- To obtain the Mbh — -Mbuigc correla- 
tion in normal galaxies, the observing sequence is to first select bulges or galaxies from 
^(■Mbuige), the galaxy bulge mass function shown schematically in Fig. [Sf>, followed by 
measuring the BH through stellar dynamics or other means. The latter probability: mea- 
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Figure 5. Measuring BHs in normal galaxies, a) Given the intrinsic prior of Fig. |U what is 
actually observed when objects are first drawn from the bulge mass (i.e. luminosity) function 
^Buigc (Panel b), followed by perfect Mbh measurement. Note that {M bulge) (with Mbh as 
the independent variable, as represented by horizontal rectangles as visual aid) trends upward 
(dashed line) even though {Mbh} (solid line, with .Mbuige as the independent variable) is not 
biased. 



Figure 6. Measuring galaxies around luminous quasars. Selecting quasars from the quasar mass 
(i.e. luminosity) function <1?bh, given that the intrinsic correlation P(A^bui gc | A^bh) between 
Mbh and M bui gc is linear (Fig.rjJ. The tapering of the correlation at high .Mbh is due to there 
being fewer luminous quasars in the universe, as illustrated by the mass function $(7Mbh) to 
the left. Note that {Mbh} (at a given .Mbuige) trends to the right (dashed line), but (Mbuige) 
(solid line) is not biased. 

suring a BH of mass Mbh after selecting on A'fbuige) is the conditional P(A / JBH|A / (buigc) 
of Fig. 01 The observational sequence: $(Afbui g e) x P(AfBH| Albuigc) therefore estab- 
lishes the observed correlation P(A / (BH,A / (buigo) shown in Fig.JSJa. The effect of selecting 
on $(A / /buige) tapers off the underlying correlation at the right side indiscriminantly of 
■Mbh- This LF imprint is present even if every galaxy and BH can be detected and 
measured precisely. It is not a Malmquist bias, and the effect prevents us f rom directly 



observ ing the intrinsic correlation. Fig.[5]a is the same, in essence, as Fig. 2 in lLauer et al 
( 20071 ). The tapering by $(A1buigc) causes the (.Mbuige) f° r a given ./Mbh to deviate from 
the intrinsic trend, as illustrated by the dashed line in Fig.[5]a. It leads to the notion that 
"at a given BH mass, there are more low mass host galaxies than high mass." This is the 
main reason behind the argument that high-z data in Fig. [3] are biased. 

However, what is subtle and widely misconstrued is the fact that the distribution of 
Fig. [Sja applies only to normal galaxies but not for quasars. In quasars, the reverse 
observational sequence occurs, i.e. measuring host galaxies around BHs, as opposed 
to measuring BHs in normal galaxies. When selecting on quasars, there is an agree- 
ment that one does not draw them from the galaxy luminosity function, but instead 
from $(Lqso) _ th e quasar luminosi ty function. Moreover, because A^bh in quasars 
scales like £qso 1 Kaspi et al. Il200fi. and quasars appear to radiate at fixed fraction 
of Eddington ratio (jKollmeier et al.l 120061) . selecting on £qso is essentially drawing on 
< I > (A / (BH)(Fig.[U left). After selecting on quasars, the host galaxy masses Albuige are then 
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drawn from the conditional probability of finding .Mbuigc around a BH of mass .Mbh, i-c. 
P (-M bulge | .Mbh)- Therefore, the observational sequence for quasars is given by the prod- 
uct: P(-M B H,A4bui g c)= ®{M B h) x P(A^bui g c|A^BH)- Comparing this to P(7W B H,A^bui g e) 
for normal galaxies, one notices the labels of .Mbh an d .Mbuigc are simply switched. Do- 
ing so leads to Fig. [13 whereby the ^(.Mbh) selection attenuates the intrinsic correlation 
of Fig. [Hon the upper (.Mbh) side. In other words, in quasars it is (.Mbh), not (.Mbuige), 
which is lower than intrinsic. Therefore, the fact that high-z data in Fig. [3] lie on the 
opposite side of the expected trend is a testament to a positive evolution in T if the BH 
mass scale is absolute. However, it is not, as discussed below, which means this effect is 
only secondary. 

It appears one reason there is widespread misconception on this issue is a tendency to 
apply the intuitive notion that there are more low mass than high mass galaxies, when 
doing so is not appropriate. In other words, after selecting a quasar from the BH mass 
function $(.Mbh), the tendency is to believe the host galaxy should be drawn from the 
galaxy mass function, P(M bulge |-Mbh)= bulge), because there are more low mass 
than high mass galaxies, as opposed to the intrinsic correlation of Fig. Doing so leads to 
the joint product P(MBH,-Mbuigc)= ( E > (.Mbuige)x ( I ) (.MBH) which heuristically produces 
a distribution given by Fig. [7] Clearly, observations do not support this because the 
j oint p roduct produces no correlation between .Mbh and .Mbuigc- Note that lLauer et al " 



t p: 

( 20071 ) did not make this particular error; their conditional probability comes from the 



linear correlation of Fig. [H not $(.Mbh)- Statistically, the only way for the conditional 
P(A(buigc|A / (BH)= ^(A^buigo) is if .Mbh and .Mbuigc are intrinsically unrelated. 

3.2. The Effects of Malmquist and Quasar-to-Host Galaxy Contrast Biases on the 
Mbh — -Mbuigc Correlation in Quasars. 

Malmquist bias is another common factor used to argue against findings that the ratio of 
-MBH/-Mbuigc is higher at high-z than now. However, Fig.[5]illustrates schematically that 
Malmquist bias only attenuates the underside of the distribution. It does not affect the 
trend at the massive end. It is also qualitatively very different from high-z observations 
of Fig. [3] because, as seen in Fig. the attenuation is uniform at a constant .Mbh; it 
does not cause the points to lie systematically to the left of the correlation line. 

Lastly, measuring host galaxies around quasars is affected by the fact that only lu- 
minous host galaxies can be detected from beneath luminous quasars. Under normal 
circumstances without gravitational lensing, host galaxies of quasars are extremely diffi- 
cult to detect when the quasar:host ratio is larger than 10 : I at a seeing of 0.1 arcsec. 
This selection bias tapers the correlation along a diagonal line illustrated schematically 
in Fig. the angle of the diagonal depends on the magnification ratio if the quasar 
sample is from gravitational lenses, as in Fig. [31 Nevertheless, the observational pressure 
is to shift the (.Mbh) and (.Mbuige) averages to the right of the intrinsic correlation. 

Comparing Fig. [8] with Fig. [3] therefore qualitatively illustrates that the finding of a 
larger A^Bn/A^buige ratio in high-z quasars is not due to known luminosity selection 
effects. Qualitatively, observational pressures greatly favor galaxy detections to the right 
of the correlation line where the quasar luminosity contrast is low. The missing objects to 
the right of the intrinsic correlation may be caused by quasar surveys that fail to classify 
low contrast, thus redder, objects as being quasar candidates . Howeve r , given that even 
redder and lower contrast systems make it into the IPeng et al.l (|2006ah quasar sample at 



z = I , this effect is judged on face value to probably not be the main culprit. 

Note that, hypothetically, it is possible for studies using other selection functions beside 
ones mentioned to distill a sample of low luminosity quasars that are then found to the 
right of the correlation. That would not necessarily contradict current conclusions using 
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Figure 7. The joint distribution P(MBH,-Mbui g e)=&{-Mbui g e)x&{MBti)- A common, but 
wrong, notion that, for every quasar observed from the BH mass function ($bh), the host galax- 
ies can thereafter be drawn from the bulge mass function $Bui g c = P(-Mbuigol-MBH)- Rather 
than producing the Mbh — .Mbuige correlation, the joint product results in no correlation. 

Figure 8. Other selection biases affecting quasar observations. Observing the Mbh — .Mbuige 
correlation in quasars has additional biases shown, given that the intrinsic correlation 
P(-Mbuigc|-MBH) between .Mbh and .Mbuige is linear (Fig. [4j . The tapering of the correlation at 
high .Mbh is again due to selection on $(.Mbh). Malmquist bias selects against faint quasars, 
independent of .Mbuige (lower dashed line). On the other hand, faint host galaxies sitting be- 
neath luminous quasars are hard to detect, giving rise to a diagonal selection bias. The exact 
angle of the diagonal bias depends on the degree of host galaxy magnification. Both (.Mbh} 
and (A4 buigo) are now shifted to the right of nominal center (solid line) due to contrast bias. 
Comparing this expected distribution with Fig.|3]shows a lack of LF bias in high-z observations. 

quasars. Instead, that hypothetical sample can have properties that distinguish them 
physically from the host galaxies of luminous quasars. Selection functions that draw on 
different physical attributes may find objects in different parameter space of the same 
underlying .Mbh — ■Mtotal correlation. This might explain the di fferent conclusions see n 
between quasars and sub-mm galaxies hosting active nuclei (e.g lAlexander et al.lf2008l) . 
To talk about evolution, it is therefore necessary to compare objects selected based on the 
same physical and observational selection functions. In that respect, the quasar-quasar 
comparisons of high-z and low-z are currently the most internally consistent sample to 
address the issue of the AIbh — A^buigc evolution. Lastly it is important to note that 
where selection biases strongly partition observable parameter spaces, it is important to 
not only consider the mean of some trend, but also the distribution as a whole. Taking 
out biases in distributions from known selection functions are feasible. 



3.3. The Black Hole Mass Scale in Quasars Is Tied to Normal Galaxies Through the 

.Mbh — o~* Correlation 
In the context of the evolution in AffiH/A'fbuige ra ti° L, the discussions above on the 
luminosity function bias is mostly academic because th e BH mass scale in quasars is tied 
to normal galaxies through the .Mbh — er* correlation (jOnken et al.ll2004l ). The bias due 
to the LF selection is normalized out to first order. 

To second order, there are other concerns when comparing the high-z sample with 
low-z, such as the relative luminosities of the quasars, the Eddington ratio, and the 
possibility that high-z BHs are unusually massive. These concerns are addressed by Fig.|9l 
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Figure 9. Comparison of quasar properties at high and low z. a) Quasar radiating efficiency 
in units of [Xqso /i/Eddington] ■ 6) .Mbh in the quasar sample in LW©]. c) Quasar luminosity in 
[ergs s ]. The solid line and dashed line are reference lines showing th e trend in the evolutio n 
of the quasar break luminosity L£oi (with arbitrary normalization) from I Hopkins et all (|2007bT l. 
Note that it is the ratio of AL = log(Lv/L^ ol ) that affects the degree of bias, so high-z quasars 
are not relatively biased despite the Lqso being somewhat higher. Note that high-z qua s ars ar e 
not atypically large in all these observables compared to low-z quasars from lKim et al.l (2008). 



which shows that the high-z sample is not too different from the low-z sample in those 
respects. The one caveat is that, even though the systematic bias in LF selection above 
is normalized out to first order, there remains residual biases relative to some reference 
point of the .Mbh — -Mbuigc correlation. Objects more, or less, luminous compared to 
that reference point may lie systematically away from the correlation, keeping in mind 
this is at most a second order effect. Taking that pivot point to be around t he break of 



the $(A 4bh), one can see in Fig. [5b. that the high-z quasar luminos ities in the Peng et al 



(|2006al) sample track the evolution of the LF break (taken from iHopkins et al.l 12007b 



with arbitrary normalization) of the quasars fairly closely both at low and high redshifts. 

The fact that the .Mbh scale in quasars is normalized to normal galaxies means that 
claims of .Mbh — -Mbuigc evolution is only meaningful if low-z quasars do not show the 
same offset. Figure [TUb shows that the low-z quasars scatter around the normal galaxy 
correlation (solid line), which indicates the bootstrapping does not leave large residual 
biases. In contrast, the high-z sample Figure [TUb clearly lies off the correlation, despite 
the .Mbh, luminosity, and Eddington ratio being quite similar to the low-z sample. 



Statistics of Mbh — M. bulge Correlation and Co-evolution 



9 



M„ (total) [mag] 




log (M, M /M ) (total) log (800 M BH /r/M Q ) 

Figure 10. (On left) Comparing the A^bh — Albuigc relation for low-z and high-z quasars hosts, 
a) Low-z quasars. There is no offset relative to normal galaxies (solid line) because the BH mass 
scale in quasars is normalized to agree, b) However, high-z quasar hosts exhibit an offset. 

Figure 11. (On right) The radius (r e ) vs. "bulge mass" relation of high-z quasar hosts compared 
with low-z elliptical galaxies. The solid circles are normal elliptical galaxies at z = with 
dynamical .Mbh meas urements. The squa re data points come from gravitationally lensed quasar 
host measurements of lPeng et alj f2006a). The product 800.Mbh /T is the expected bulge mass 
at the observed epoch as inferred from Fig. [3] or the right hand side of Fig. 1101 Quasar hosts at 
z > 2 are very compact. 



4. Quasar Host Galaxies at z > 2 Are Under-Sized for Their Mass 

Another interesting evidence that high-z quasar hosts have a larger T =Mbh/M bulge 
relative to z — (which can also be thought of as a mass deficit in the bulge) comes from 
comparing the size vs. A^buigc correlation at the observed epoch with galaxies today, as 
shown in Fig. 1111 In that Figure, the host galaxy mass is inferred from the luminosity 
of the host galaxy. But, it is useful to recast the mass in terms of T, so as to emphasize 
how the controversial mass deficit parameter affects the size- M bulge correlation in high- 
z quasars. Doing so, the host galaxy bulge mass is: A^buige(^) = 800x Mbh/F. This 
equation comes from the fact that normal galaxy bulges at z = have T = 1 and 
A^buigcfz = 0) = 800xA4bh- The correlation of r e with A/fbuige is revealing because 
unknown luminosity selection biases are effectively normalized away by accounting for T. 
Figure [TT1 shows that the host galaxies at z « 1 seem to lie on the size-A-^buige correlation, 
whereas higher redshift host galaxies appear to be much more compact per unit mass. 
By z « 2, the host galaxies appea r to be too small by a factor of 2-3 compared to normal 
galaxies of the same mass today riPend 120041) . One way to weaken the conclusions is for 
r to be even larger than the controversial claim, which permits these objects lie on the 
modern day size-mass correlation. The fact that massive galaxies at high z appear to b e 
very compact has been ob s erved by a num b er of studies, including Truiillo et alj ( 20061 ); 



van Dokkum et ail (|2008l ); IStockton et al. ( 2008 ) among others, and may point to the 



same evolutionary pathways between quasar host galaxies and distant red galaxies. 
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5. Conclusion 

The statistics of galaxy merging may shed some light on the co ntroversial fi nding of an 
evolution in the A4bh/ Alga] ratio. In Monte-Carlo simulations of lPend ( 2007 ). high mass 



objects often tend to lie to the left of the asymptotic linear correlation. This happens 
because such objects were large outliers in the initial distribution, thereby taking them 
longer time to evolve onto the asymptotic relation. Another potential explanation for 
the larger Mun/Mgai ratio is that the quasar phase may signify recent BH growth, so 
by observing luminous quasars we catch them in a special s tate on the A^bh — A fhuip-R 



correlation. This is consistent with the simulation findings of lHopkins et al. ( 2007a ) who 



explain large offsets as being due to gas rich mergers that both feed the central BH and 
possess a larger mass fraction in gas. As explained by merger statistics, the temporary 
up-tick in the BH mass can subsequently merge back onto the asymptotic linear correla- 
tion through minor mergers. Indeed, this is seen in Monte Carlo simulations where the 
BHs were artificially boosted in mass followed by regular mergers. Combining statistical 
simulations with observations that high-z quasar host galaxies are very compact, and 
the fact that major mergers do not change the A'fBH/A'fbuige ratio, seem to consistently 
point to minor mergers being important for transforming quasar hosts morphologically 
from their compact state at z ~ 2 into massive, extended, elliptical galaxies today. 

I thank Jenny Greene, Chris Kochanek, and Luis Ho for providing very thoughtful com- 
ments. 



References 

Adelberger, K. L., & Steidel, C. C. 2005, ApJL, 627, LI 

Alexander, D. M., Brandt, W. N., Smail, I., & Swinbank, A. M., et al. 2008, AJ, 135, 1968 

Ciotti, L., & van Albada, T. S. 2001, ApJL, 552, L13 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604 

Ferrarese, L., & Merritt, D. 2000, ApJL, 539, L9 

Fine, S., Croom, S. M., Miller, L., Babic, A. et al. 2006, MNRAS, 373, 613 
Gebhardt, K., Bender, R., Bower, G., Dressier, et al. 2000, ApJL, 539, L13 
Granato, G. L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 2004, ApJ, 600, 580 
Haxing, N., & Rix, H.-W. 2004, ApJL, 604, L89 

Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., & Krause, E. 2007a, ApJ, 669, 67 
Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007b, ApJ, 654, 731 
Islam, R. R., Taylor, J. E., & Silk, J. 2003, MNRAS, 340, 647 

Kaspi, S., Smith, P. S., Netzer, H, Maoz, D., Jannuzi, B. T., & Giveon, U. 2000, ApJ, 533, 631 

Kauffmann, G., & Haehnelt, M. 2000, MNRAS, 311, 576 

Kim, M., Ho, L. C, Peng, C. Y., Barth, A. J. et al. 2008, ApJ, 687, 767 

Kollmeier, J. A., Onken, C. A., Kochanek, C. S., & Gould, A., et al. 2006, ApJ, 648, 128 

Kormendy, J., & Richstone, D. 1995, ARAA, 33, 581 

Lauer, T. R., Tremaine, S., Richstone, D., & Faber, S. M. 2007, ApJ, 670, 249 
Magorrian, J., Tremaine, S., Richstone, D., Bender, R. et al. 1998, AJ, 115, 2285 
Marconi, A., & Hunt, L. K. 2003, ApJL, 589, L21 

McLure, R. J., Jarvis, M. J., Targett, T. A., Dunlop, J. S., et al. 2006, MNRAS, 368, 1395 
Onken, C. A., Ferrarese, L., Merritt, D., Peterson, B. M., et al. 2004, ApJ, 615, 645 
Peng, C. Y. 2004, PhD thesis, The University of Arizona, United States - Arizona 
— . 2007, ApJ, 671, 1098 

Peng, C. Y., Impey, C. D., Ho, L. C, Barton, E. J., & Rix, H.-W. 2006a, ApJ, 640, 114 
Peng, C. Y., Impey, C. D., Rix, H.-W., Kochanek, C. S., et al. 2006b, ApJ, 649, 616 
Robertson, B., Hernquist, L., Cox, T. J., Di Matteo, et al. 2006, ApJ, 641, 90 
Salviander, S., Shields, G. A., Gebhardt, K., et al. 2006, New Astronomy Review, 50, 803 



Statistics of A4bh — A'lbuigo Correlation and Co-evolution 

Stockton, A., McGrath, E., Canalizo, G., Iye, M., & Maihara, T. 2008, ApJ, 672, 146 
Treu, T., Woo, J.-H., Malkan, M. A., & Blandford, R. D. 2007, ArXiv e-prints, 706 
Trujillo, I., Forster Schreiber, N. M., Rudnick, G., Barden, M. et al. 2006, ApJ, 650, 18 
van Dokkum, P. G., Franx, M., Kriek, M., Holden, B., et al. 2008, ApJL, 677, L5 
Woo, J.-H., Treu, T., Malkan, M. A., & Blandford, R. D. 2006, ApJ, 645, 900 



